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O ■ ABSTRACT 

^ r"! , We present the first detailed structure formation and radiative transfer simulations of the reion- 

Q^. ization history of our cosmic neighbourhood. To this end, we follow the formation of the Lo- 

cal Group of galaxies and nearby clusters by means of constrained simulations, which use 
the available observational constraints to construct a representation of those structures which 
reproduces their actual positions and properties at the present time. We find that the reion- 
^ I ization history of the Local Group is strongly dependent on the assumed photon production 

efficiencies of the ionizing sources, which are still poorly constrained. If sources are relatively 
efficient, i.e. the process is 'photon-rich', the Local Group is primarily ionized externally by 
the nearby clusters. Alternatively, if the sources are inefficient, i.e. reionization is 'photon- 
poor' the Local Group evolves largely isolated and reionizes itself. The mode of reionization, 
external vs. internal, has important implications for the evolution of our neighbourhood, in 
terms of e.g. its satellite galaxy populations and primordial stellar populations. This therefore 
^vf^ . provides an important avenue for understanding the young universe by detailed studies of our 

• ■ nearby structures. 

^^ I Key words: H II regions — ISM: bubbles — ISM: galaxies: halos — galaxies: high-redshift — 

^^ . galaxies: formation — intergalactic medium — cosmology: theory — radiative transfer — meth- 

ods: numerical 

•i-H 

^ ■ 1 INTRODUCTION This strongly suggests that the reionization epoch started well be- 

' ' ' , . . , . , .„. , . . , , fore redshift 10 and therefore was fairly extended in time. It also 

Approximately thirteen billion years ago the cosmic neighbour- ^ , . , j .1 .1 • . r • • .• 1 

, ^, , . / „ , -^ * , ^ r' , confirmed independently the existence 01 a reionization epoch, re- 

hood destined eventually to become our Local Group of galax- • j i_ .1 • j j-.- 1 .• 1 j .1 ry., r •• .• 

/ . . . . . . ^ quired by this additional optical depth. The process 01 reionization 

les underwent a dramatic transition: a giant ionization iront swept ,jr i- ri .ir.- 

,^ . . ^. . ,. . ^, . , had lar-reachmg consequences lor subsequent galaxy lormation. 

through, engulfing it m a sea 01 ionizing radiation. This occurred rr., 1 . • • .• 1 .• 11 • •• .• 

i\r- r It, •• r,- , • Tfic photoiomzation heating which accompanies reionization m- 

as a local manifestation of a global transition of the intergalactic , ,, , . r .1 1 r j r t^ 

., ,, . ^ ^ , ^ . .. creased the gas temperature irom the very low one, 01 order a lew K 

medium m the whole universe referred to as Cosmic Reionization, , j-.i^ t^ia ir.ix^^.r j 

, ^ , ,. . ^ , ^ , . ^, . or less, during the Cosmic Dark Ages before the first stars lormed, 

caused by the radiation irom the first galaxies. This process con- , ir^4 t^ t., • • , • , ,, j- t 

, , . , . ,. ^ ^ , 1 ,1 1 • to ^ 10 K or more. This m turn increased the corresponding Jeans 

verted the intergalactic medium from neutral and cold gas during ^, u u- u ^ .en 

, ^ -T^T* , r .• t-, • mass, the mass above which gas pressure cannot sucessiully coun- 

the Cosmic Dark Ages beiore any galaxies existed, into a hot, ion- ^ ^ . , 1 . r j r •. j ^^i • . 1 

. ^ -^ ° teract gravity, by about 5 orders 01 magnitude. This strongly sup- 

^ ' pressed the formation of low-mass galaxies and cut off the star for- 

The absorption spectra of OSOs from redshift to about 6 ^. . • 1 r j j ^u u u u u • .c 

\ \ , t^ , 1 , P „ • mation m previously-iormed ones, and thereby should have signii- 

show that the intergalactic medium has been almost fully ion- . ^1 • r, t ^1 1 1 ^- n ^ ^ ^ nv 1 • 

^ , ,.^ . ^ , ^ . ^ , , , , icantly influenced the early population 01 dwari satellite galaxies, 

ized tor most 01 the liietime 01 the Universe. On the other hand, t- .1 • •• .• • r. • 1 j 1 ui 1 

For this reason, reionization is oiten invoked as a plausible expla- 
the recent data from the Wilkinson Microwave Amsotropy Probe ^ .c ^u u j 1 1 r 1 ^ nv j ^ ^u 

. , , , , , • , 1 , r nation tor the observed lack 01 galaxy satellites compared to the 

(WMAP) satellite yielded a rather large optical depth for scatter- , j- ^ j u j 1 ^^ • 1 ^- /t^ n 1 TTI 

: , ^ . ,,. ^ , t , n . numbers predicted by pure dark matter simulations (Bullock e t al.l 

mg the Cosmic Microwave Background photons on tree electrons. U^nnl k/r ~ mUr^kn' TT ^^ — TTr^^^TTTT? "7! — \ . 1 o^TTvT 

^ & F 120001 : iMunoz et al.l 12009: Maccio et al. 2010: Busha et al. 2010). 

The same ionization and heating process also changed significantly 

the character of star formation, as stars whose formation starts out 

^ e-mail: LTIliev@ Sussex ac.uk from hot, ionized gas, even those with a primordial element abun- 
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danc e, are different from the ones forming out of initially cold gas 
(e.g. lMackevetalJl2003h . 

In recent years a variety of theoretical and numerical mod- 
elling work has shown that reionization was highly patchy in na- 
ture, with a very large variatio ns in the times at which different 
regio ns become reionized (e.g. Illiev et al .1120061 . 120071 : IZahn et al.l 
l2007h . Many of the effects of reionization on the later structure 
formation are highly dependent on the stage at which the forming 
structures are at the time when they become reionized. Therefore, 
it is important to follow the reionization history of different patches 
in detail and to compare the timing of reionization with the stage 
of f ormation of the local stru ctures. 

IWeinmann et al.l (120071) used large-scale reionization simula- 
tions to study the reionization history of galaxies of a variety of 
present-day types: field galaxies, cD central cluster galaxies, L* 
galaxies and Local Group-like systems. In particular, they focused 
on the question which galaxies are reionized internally (i.e. by 
their own progenitors) vs. externally (i.e. before sigificant fraction 
of their mass is in collapsed objects). For most galaxy types the 
answer is statistical, with certain probabilities for each outcome. 
They found that there is a halo mass scale, of order 10^^ M0, 
which divides the massive galaxies which are predominantly inter- 
nally reionized and lower-mass on es for which the opposite is true. 
More recently, I Alvarez et al.l (12009.) studied the same problem us- 
ing larger-scale, coarsely-resolved semi- analytical reionization cal- 
culations. They achieved better statistics for the larger halos, from 
Milky Way- sized up to galaxy clusters, as a consequence of the 
larger volume their calculation follow ed and reached very similar 
conclusions to lWeinmann et al.l (l2007h . 

However, these previous studies are statistical and therefore 
can only yield a certain probability for a given type of system to be 
externally or internally reionized. In contrast, the answer for a spe- 
cific system, like our own Local Group of galaxies, consisting of 
the Milky Way, Andromeda, M33 and their satellite galaxies, will 
depend on the details of its and its neighbouring systems' forma- 
tion, the timing of that formation, and the relative positions in space 
of their projenitors at the relevant epochs. This kind of detailed in- 
formation can only be obtained through numerical simulations with 
constrained initial conditions. 

Constrained simulations aim to reproduce the spatial and ve- 
locity structure of our Local Group and its neighbourhood at 
present (redshift z — 0). In this work we combine the most ad- 
vanced constrained realizations of the Local Group available at 
present with an accurate treatment of the radiative transfer during 
cosmic reionization. This allowed us for the first time to calculate 
the specific reionization history of all Local Group progenitors, as 
well as those of the nearby structures like the Virgo and Fornax 
clusters of galaxies. 



2 METHODOLOGY 

2.1 Constrained simulations of the local universe 

The optimal way of constructing a numerical simulation that 
closely repr oduces our local cosmol ogical neighborhood is pro- 
vided by the taoffman & Riba^ (Il99lh algorithm for making a con- 
strained realizations of Gaussian random field. Given our ability 
to extract observational data that can be imposed as linear con- 
straints on the primordial perturbation field, this method can be 
used to co nstruct initial conditi ons that obey these constraints. Here 
we follow iKlvpin et al.l (120031) and impose two types of data on 
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Figure 1. Collapsed fraction in high-mass (M > 10^ Mq; red, solid) and 
low-mass sources (M < 10^ Mq; blue, dashed) (insets show the same in 
log scale) vs. cosmic redshift for WMAP5 (left) and WMAP3 (right) cases. 

the simulation. The first data set con sists of peculiar vel ocities of 
galaxies, drawn from t he MARK III dWillick et al. I [19971) . surface 
brightness fluctuations ("Tonrv et al.^'2001') catalogues and the Cat- 
alog of Nearby Galaxies ( Karachentsev et al. 2004). The other data 
set is obtained from the catalog of nearby X-ray selected clusters 
of galaxies (Reiprich & Bohringer 2002). The main obstacle faced 
here is how to translate the present epoch observables into quanti- 
ties that are linear in the primordial perturbation fields. Peculiar 
velocities evolve more slowly than the density field and are as- 
sumed here to be linear. The present epoch virial parameters of 
the clusters are processed by the spherical top-hat model to pro- 
duce their linear overdensity. The data used here and its associated 
observational errors effectively constrain the LSS on scales larger 
than ^ 5h~^Mpc (cf. Klypin et al. 2003). The simulation used 
here is designed to reproduce the Local Supercluster, harbouring 
a Virgo-like cluster. Such a configuration is easily reproduced by 
the constrained simulations. Local Group (LG) like objects, on the 
other hand, are randomly emerging in the simulations. The simula- 
tions used here were each selected out of a few to have a LG-like 
object, similar to the observed one. They are described in more de- 
tail in lGottloeber et all (l2010h . 

Starting fron i the above constrained initial conditions, the 
GADGET-2 code ('Springei'2005) was used to follow the dark mat- 
ter in a L = 64h~^Mpc computational box, spanned by 1024^ 
particles, starting at redshift z = 100. The cosmological parame- 
ters given by the WMAP 3 -year data have been adopted (Qm = 
0.24, Qa = 0.76, h = 0.73, Vtb = 0.0418, ag = 0.75, n = 0.95), 
giving a particle mass of moM = 1.63 x 10jh~^M c:). (A more 
detailed description is provided in lZavala et al.l (l2009l) ). 

In order to check the robustness of our results with respect to 
the particular realization of the constrained simulation, we also per- 
formed a second simulation with an independent underlying ran- 
dom realization. This second simulation implements the same con- 
straints on the local structures at the present time. The simulation 
volume and number of particles are the same as above, while the 
background cosmology is now based on the WMAP 5-year data, 
combined with the available constraints from the large-scale struc- 
ture (BAO) and supernovae (Qm = 0.279, Qa = 0.721, /i = 
0.73, Ob = 0.046, as = 0.817, n = 0.96). In this case we in- 
clude in the analysis not just the Local Group and Virgo, but also 
our other nearby galaxy cluster, Fornax. In Figure [2 we show the 
local density distribution at redshift z — 9 within a slice of comov- 
ing 30/i"^ Mpc X 30/i"^ Mpc and 7h~^ Mpc depth which con- 
tains the projenitors of the objects of interest, Virgo, Fornax, M31 
and Milky Way. The slice is situated in the supergalactic YZ-plane. 
The main progenitors of these objects have at redshift z = 9 masses 
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Figure 2. Matter distribution in a volume of 30/i~^ Mpc x 30/i~^ Mpc x 
7h~^ Mpc depth in comoving Mpc at redshift z = 9. The circles with 
decreasing radius according to the decreasing mass of the objects denote the 
regions where the progenitors of Virgo, Fornax, M3 1 and MW are situated. 
The slice is situated in the supergalactic YZ-plane. 



of 1.2 X lO^^/i"^ Mo, 6.6 X 10^°/i"^ M©, 3.8 x lO^/i"^ M©, 
and 3.4 x 10^ /i~^ M©, respectively. The marked Virgo (Fornax) 
regions are populated with about 2000 (700) cluster progenitors, 
each with a mass larger than 4 x l(fh~^ Mq. This leads to a total 
progenitor mass of 2.4 x 10^^/i~^ M© (0.9 x 10^^/i~^ M©). One 
can clearly see the low density regions between the progenitors of 
the clusters and the progenitors of the Local Group. There are a 
few other objects next to the Local Group but with much smaller 
masses than that concentrated in the proto-cluster regions. 

All our simulations were performed within the CLUES (Con- 
strained Local UniversE Simulations) projec{j. 



2.2 Radiative transfer simulations 

The radiative transfer simulations were performed using a ra- 
diativ e transfer and non-equ ilibrium chemistry code called C^- 
Rav (iMellema et aL l2006al). tested in d etail as d iscussed in 
iMellema et aP (l2006ah . lniev & et all (l2006h an dllliev e t al. (2009). 
Our simulation methodology was presented in llliev e t al. (2006); 
IMellema et aP ([2006b) and Iliev et all ([2007h . The underlying N- 
body structure formation simulations discussed in the previous sec- 
tion provide a time sequence of density distributions and catalogues 
of identified halos. Gas is assumed to follow the dark matter distri- 
bution, which at these scales (0.25 /i~^Mpc— 64 /i~^Mpc) is a very 
good approximation. We have produced 53 density slices and halo 
catalogues roughly equally spaced in time between the redshifts 20 
and 6, every ^ 14.5 Myr. All halos are assumed to host galaxies 
and thus to be potential sources of ionizing radiation. Each is as- 
signed an ionizing luminosity proportional to its mass. Low-mass 
sources, with total mass below 10^ M©, are assumed to be active 
in the neutral regions, but to be suppressed in the ionized regions. 



Therefore, once their cell is ionized their emissivity is set to zero. 
These assumptions are based on the fact that the Jeans mass rises 
significantly when a region is ionized and heated, which thereby 
limits the fresh gas infall onto halos and supress es the future forma- 
tion of low-mass galaxies, see llliev et al.! (120070 for a more detailed 
discussion. 

Both the density and halos are binned on a 256^ grid for the 
radiative transfer processing. Halos in the same RT cell are com- 
bined and a luminosity is assigned following a simple, physically- 
motivated prescription. The emissivity of the ionizing sources - 
how many ionizing photons produced in galaxies over some time 
period reach the intergalactic medium - depends on what fraction 
of the galactic gas is converted into stars, how effective are the stars 
at producing ionizing photons and, finally, what fraction of the pho- 
tons manage to escape the galaxy. We parametrise it with a single 
parameter /^, which is equal to the number of ionizing photons 
per atom in galaxies which reach the intergalactic medium in the 
time between two consecutive time slices The main difference be- 
tween our two adopted background cosmologies is in the density 
fluctuation normalization (as). A higher as yields higher halo col- 
lapsed fractions at any given epoch and a correspondingly larger 
number of collapsed halos, as illustrated in Figure [T] In terms of 
reionization this means that for a fixed ionizing emissivity propor- 
tional to the collapsed fraction the evolution is shifte d to somewhat 
earlier times, resulting in an earlier overlap epoch (lAlvarez et alj 
12006). This effect was much larger for WMAP 1-year (as = 0.9) 
vs. WMAP 3-year data (as = 0.74) than it is for WMAP 5-year 
data ((78 = 0.8) vs. the WMAP 3-year data. 

Much larger uncertainty concerns the properties of high red- 
shift sources. Because of the still scarse observational data those 
are not at well constrained at present and hence /^ is largely a 
free parameter within certain, fairly wide, bounds. Furthermore, the 
radiative transfer simulations are quite computationally expensive, 
and therefore it is not practical to investigate the full available pa- 
rameter space. We address this problem by considering two repre- 
sentative cases which are roughly bracketing the range of expected 
behaviour, while at the same time satisfying the available global 
observational constraints, as follows. For our Model 1 simulation 
we use the first N-body realization, with WMAP3 cosmology and 
we adopt /^ = 100 for the massive sources (M > 10^ Mq) and 
/^ = 250 for the low-mass, suppressible ones (M < 10^ Mq). 
These relatively high efficiencies yield a fairly fast, photon-rich 
reionization process with an early overlap. For our second simula- 
tion, on the other hand, we use the second constrained realization, 
with WMAP5 cosmology, and we adopt a lower source efficien- 
cies of /^ = 10 for the massive sources and f^ = 150 for the 
low-mass ones. This yields a more extended, photon-poor reioniza- 
tion history. For conciseness of notation, we will to refer to the two 
combinations of N-body realizations and corresponding radiative 
transfer simulation simply as the Model 1 and Model 2 cases, with 
the implicit understanding that this implies a different constrained 
realization and assumed photon emissivities of the sources, as well 
as different background cosmology. 



^ One can also introduce a slightly different efficiency parameter, g^, given 
by 5^7 ^ 



/^ i — ^^ j where At is the time between two snapshots from 
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the N-body simulation. This has the advantage that it is a rate per unit time 
and as such it is independent of At, which makes easier comparisons be- 
tween simulations with different At. 
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Figure 3. (Bottom) Evolution of the mass-weighted ionized fractions, Xm 
for Model 1 (red, dashed) and Model 2 (blue, solid) cases, inset shows the 
same in linear scale; and (top) the corresponding ratios of mass-weighted 
to volume- weighted ionized fractions, which corresponds to the average 
density of the ionized regions in units of the mean, versus redshift z. 



3 RESULTS 

3.1 Global reionization history 

The global mass-weighted reionization histories produced by our 
two simulations are shown in Figure [3] In either case reionization 
starts when the first resolved halos (which correspond to the first 
ionizing sources) form, around z ~ 20. 

Our adopted /^ values yield final overlap of ionized regions 
at Zov ^ 9 (8) for the Model 1 (Model 2) case, in rough agree- 
ment with the current observational constraints. The correspond- 
ing integrated Thomson scattering optical depth seen by the Cos- 
mic Microwave Background photons is res = 0.094(0.069) is 
also in agreement with the latest constraints from WMAP satellite 
combined with the other available datasets, Tes = 0.084 ± 0.016 
(iKomatsu et al.''2009). 

The early reionization (z > 14) is driven primarily by the low- 
mass sources, which have similar efficiencies in the two cases (the 
slightly lower source efficiency in the Model 2 case is compensated 
for by its higher collapsed fraction at the same redshift) and as a 
consequence the two reionization histories are initially very similar. 
Later on the larger sources take over, both because of their rapidly 
rising collapsed fraction (cf . Figure[T]) and the strong Jeans suppres- 
sion of the low-mass sources, and thus reionization proceeds more 
slowly in the Model 2 case due to the lower efficiency adopted for 
its high-mass sources. The fact that the mean mass-weighted ion- 
ized fraction, Xm is always larger than the corresponding volume- 
weighted one, Xv (Figure[3] upper panel) indicates that reionization 
proceeds in an inside-out fashion (i.e. high-density regions are pref- 
erentially ionized first) in both cases, in agreement wit h previous 
sim ulation results based on non-constrained realizations jlliev et al.l 
I2OO6) . For the Model 2 realization this ratio is noticeably higher, up 
to ~ 3 at 2; = 20, due to its more advanced structure formation at 
any given redshift. This also boosts the dumpiness of the gas, and 



therefore the recombinations, which in turn extends reionization 
even further. 



3.2 Local reionization histories 

The mean reionization history presented above ensures that the 
currently available global observational constraints - the electron- 
scattering optical depth and overlap epoch are satisfied. However, 
to achieve our present goals we need to track separately the reion- 
ization history of the progenitors of each object of interest, namely 
the Local Group, as well as the nearby clusters of galaxies. To this 
purpose, we extracted the Lagrangian mass distribution for each 
object (i.e. the mass which eventually will end up in that object by 
the present day) and followed the reionization history of all radia- 
tive transfer cells containing at least one particle which ends up in 
that object by z = 0. The resulting local reionization histories are 
shown in FigurelD along with the global one for direct comparison. 

The Lagranian regions of both the LG and the nearby clus- 
ters are significantly overdense in either constrained realization and 
at all times, reflecting the fact that all of these objects correspond 
to high peaks of the density field. The proto-Local Group region 
starts only moderately overdense, by ^ 7% (^ 9%) in the Model 1 
(Model 2) case, which rises over time as the corresponding object 
collapses gravitationally, to reach ^ 16% (^ 25%) by the global 
overlap epoch at z ^ 9 (8). The proto-clusters correspond to still 
higher peaks of the density field. Initially the proto- Virgo region is 
overdense by 12% for both simulations, rising over time to 24% 
(30%) for Model 1 (Model 2). The proto-Fornax region (Model 
2) starts 10% overdense, rising to 26% by the global overlap. The 
higher local density yields an (exponentially) larger halo collapsed 
fraction and thus ionizing photon production. Therefore, for both 
LG and clusters we can expect local reionization to occur earlier 
than average, which is confirmed by our simulation results (Fig- 
ure |4]). 

In both simulations the LG reionization starts at about z ^ 
12.5, at which time its oldest progenitor halos form. Before z ^ 
12.5 the LG ionized fraction is tiny, below 3xl0"^(2xl0"^) 
for Model 1 (Model 2). Thereafter the (proto-)LG reionization ac- 
celerates, albeit only gradually. For Model 1 (Figure ID left) the LG 
ionized fraction reaches \7% by z ^ 12 and 39%byz ^ 10.5. 
After that point the evolution becomes very fast and full ionization 
(xm > 99%) is achieved by z ^ 10. In contrast, the reionization 
history of the (proto-)Virgo cluster in the same Model 1 simula- 
tion is quite different. Proto- Virgo is a higher density peak and the 
formation of the local nonlinear structures is therefore accelerated. 
Hence the local reionization proceeds faster, as well. The evolution 
remains smooth throughout, with no sudden changes of slope, un- 
like in the LG case. The mass-weighted ionized fraction reaches 
ll%byz = 12.6, 50% by z = 11.3, and 96% by z = 10.5. 

The reionization histories are similar in the Model 2 case (Fig- 
ure ID right). Once again, the proto-cluster regions, both Virgo and 
Fornax, reionize earlier than LG and much ealier than an average 
region - Xm = 0.1 is reached by 2; = 12.3 (12.6), Xm = 0.5 
by z ^ 10.25 (10.5) and Xm = 0.9 by z ^ 9.2 (9.3) for Virgo 
(Fornax). For both proto-clusters overlap (xm = 0.99) is reached 
at 2; = 8.7 and the evolution remains smooth throughout. Interest- 
ingly, most reionization stages (but not the local overlap, which is 
roughly simultaneous) of the Fornax reionization occur earlier than 
the corresponding ones for Virgo, even though Fornax has lower 
mass at the present epoch. In comparison, the reionization of the 
LG occurs later, reaching Xm = 0.1 by 2; = 11.6, Xm = 0.5 by 
z = 9.4, Xm = 0.9 by 2; = 8.6 and local overlap is achieved by 
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Figure 4. Mass -weighted mean ionized fractions, Xm, for the Local Group, nearby clusters and global mean (as indicated by color) vs. redshift for Model 1 
(left; left to right: global, LG, Virgo) and Model 2 (right; left to right: global, LG, Virgo, Fornax). 





Figure 5. Volume renderig of the total mass, which by the present time {z = 0) will end up as part of the objects of interest here: (left) Virgo cluster (top 
object) and the Local Group (bottom object) at redshift 2; = 10; (right) Virgo cluster (top object), the Local Group (middle, smaller object) and Fornax (bottom 
object) at redshift z = 9. The color scale units are Solar masses per radiative transfer cell. 



z — ^A. It lags the global mean in its earliest stages {z > 12.5), but 
as more progenitor halos form it catches up and then speeds ahead 
after z = 12. The Local Group reionization history is again much 
less smooth than the proto-cluster ones, with significant changes of 
slope around z ^ 12, z ^ 11 and z ^ 9.6. Evolution becomes 
extremely fast after z ^ 9.6, whereby the ionized mass fraction 
jumps from 0.4 to 1 over redshift interval of just Az ^ 1. Com- 
pared to the Model 1 case the reionization histories for Model 2, 
both mean and local are ones, are much more extended in time due 
to the lower source emissivities assumed. 

Moving on to a more visual representation of the local reion- 
ization history, in Figure |5] we show a 3D volume rendering of the 
distribution of mass which is destined by the present time (z = 0) 
to become part of the Local Group and nearby clusters for Model 1 
at z = 10 (left) and for Model 2 at 2; = 9 (right). In the realization 



corresponding to Model 1 the proto- Virgo is a quite large, elongated 
object, extending for about 10 x 20 comoving Mpc, which corre- 
sponds to about 1 proper Mpc at this early time, while the proto- 
Local Group is a much smaller object a few Mpc in size, about 7 
comoving Mpc (less than 1 physical Mpc) away from proto- Virgo. 
In the alternative. Model 2 realization the size of the proto-Local 
Group and its distance from Virgo are similar to the Model 1 real- 
ization. However, the proto- Virgo is more compact, less extended 
object and a proto-Fornax is identified, as well, at a similar dis- 
tance from LG, positioned in almost diametrically opposite direc- 
tion. Proto-Fornax is a somewhat less extended object than proto- 
Virgo. 

In Figures [6] and |7] we show the evolution of the neutral mass 
remaining in the objects of interest in a time sequence for simula- 
tions Model 1 and Model 2, respectively. In each case the sequence 
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Figure 6. Evolution of the neutral mass at (top to bottom and left to right) redshifts , 
ionized. 



10.75, 10.5, 10.25 and 10 for Model 1. Red is neutral, green is 



covers the key time interval during which most of the LG mate- 
rial is ionized. In Model 1 at 2; = 10.75 the reionization of Virgo 
is already well- advanced, with a local ionized fraction above 83%, 
while LG remains largely neutral, at less than 35%. Shortly there- 
after, an ionization front arrives at the proto-LG position from the 
direction of Virgo and quickly sweeps through it as illustrated by 
the next three images. This corresponds to the dramatic jump of the 
LG local ionized fraction seen in Figure|4]at z = 10.5 — 10.25. By 
z = 10 the Local Group material is almost fully ionized, reaching 
97.5% ionized fraction by mass. In this case, therefore the Local 
Group is largely ionized from the outside, primarily by the Virgo 
progenitors. 

The evolution proceeds quite differently in the photon-poor. 
Model 2, case, as illustrated in Figure |7] While the nearby clusters, 
in this case both Virgo and Fornax, once again reionize themselves 
from the inside and relatively earlier than the Local Group, there are 
no clear ionization fronts to arrive from them and sweep over the 
LG material. The Local Group internal sources carve ionized bub- 
bles from the inside and eventually manage to reionize all the LG 
material mostly by themselves. While we cannot exclude modest 
contributions from Virgo and Fornax, the local LG sources appear 
to dominate the evolution in this case. 

In order to evaulate the local reionization process for each 
structure in a more quantitative way, we counted and added up all 
ionizing photons emitted by sources within the same Lagrangean 
volume, normalized by the total number of atoms belonging to that 



object. We also counted and added together the cumulative number 
of ionizing photons used up to ionize each object of interest and 
also to keep it ionized (i.e. recombinations, since after each recom- 
bination back to neutral state that atom would need to be ionized 
again), again normalized per atom in that object. Comparing the 
values of these two numbers over time shows if that particular ob- 
ject by itself produced enough photons up to that point in time to 
fully account for its current ionization state. Results for both Mod- 
els are shown in Figure [8] 

The results confirm our conclusions based on the visual exam- 
ination of the reionization process we discussed above. The galaxy 
clusters initially produce most of the photons needed for their own 
reionization. Until z -- 10.5 (^ -- 9 - 9.5) for Model 1 (Model 2) 
there is some deficit of photons, i.e. a little more are used up than 
produced by the progenitors of that cluster. The reason for this is 
that clusters are located at high peaks of the density field, and it 
is well established that CDM halos, and thus our ionizing sources, 
strongly cluster around such high peaks. As a consequence of that 
there are many nearby sources which surround the proto-cluster re- 
gion and contribute to its reionization. However, we note that the 
final overlap in each case is only reached after the cluster itself has 
produced sufficient number of photons. We also note that recom- 
binations have a very significant effect for proto-clusters, yielding 
usage of 1-1.5 additional ionizing photons per atom in addition to 
the one photon needed to achieve the initial ionization of that atom. 

On the other hand, the proto-Local Group results are different 
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Figure 7. Evolution of the neutral mass at (top to bottom) redshifts z 
the Model 2 model. Red is neutral, green/blue is ionized. 



9, 8.9, 8.7, 8.55 and 8.4 in 3D volume rendering (left) and cross-section (right) for 
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Figure 8. Cumulative number of ionizing photons for ionizations and recombinations (solid) and photons emitted (dashed), both per atom, for the Local Group 
and nearby clusters (as indicated by color) vs. redshift for Model 1 (left panel) and Model 2 (right panel). 



from the cluster results in both cases. In the photon-rich Model 1 
up to z ^ 10.5 there are exactly as many photons produced by 
the LG progenitors as are used up for its own reionization. After 
that point however there is a significant change of slope for both 
curves. Suddenly there are clearly more photons arriving than are 
produced locally. This change of slope occurs exactly when the ion- 
ization front from Virgo was seen above to sweep through the LG 
material. The number of photons used for ionizations and recom- 
binations is above the number of locally-emitted photons and it re- 
mains so up to the redshift of full reionization (i.e. local overlap) 
of the Local group, z ^ 10. This behaviour is more clearly seen 
in the top panel of Figures] (left), where we show the ratio of pho- 
tons used over photons emitted locally for LG and Virgo. This ratio 
starts above one at^ ^ 12.3 — 11.7, indicating that LG the early 
LG reionization is partially helped along by nearby sources that do 
not belong to it. However, this photon ratio then falls slightly below 
one and remains so until z — 10.5. At this point the external ion- 
ization front arrives and the ratio jumps to well above one, peaking 
at 1.2 and remaining above 1 until the local overlap. We therefore 
conclude that in this case most of the Local Group gas (up to 70%, 
see Figure lU was indeed reionized externally, predominantly by 
the nearby Virgo cluster. 

In our Model 2 (photon-poor) case the photon budgets of the 
proto-clusters Virgo and Fornax are very similar to each other as 
well as to that of proto- Virgo in our Model 1 . There is a noticeable 
contribution of photons from nearby clustered sources, which sup- 
ply at least 40% more photons than the respective internal sources. 
However, after z ^ 10 (when Virgo is 60% ionized and Fornax is 
70% ionized) the ratio of consumed over produced photons drops 
quickly and by their full ionization (z = 8.7) both clusters have 
internally produced more photons than are needed for their own 
ionization. The local recombinations have a significant effect, re- 
sulting in about one additional photon per atom needed to keep 
them ionized. 

In contast, the ionizing photons which are produced locally by 
the Local Group progenitors roughly balance the ones consumed 



for its own reionization all the way up to its neighbourhood overlap 
at 2; '-^ 8.4. Once again initially (at z > 11) there is some external 
contribution of ionizing photons from nearby sources, but in this 
case throughout most of the evolution (8.4 < z < 11.5) the emit- 
ted and consumed photons balance within less than 10%. In fact, at 
the later stages of the evolution the LG progenitor halos produce a 
bit more photons than are actually used to ionize it and keep it ion- 
ized. Hence, the Local Group evolves largely isolated in this case. 
After the local overlap the LG produces even more photons than are 
actually needed to keep it ionized. The effect of recombinations on 
the LG reionization history is more modest than in the case of the 
proto-clusters, resulting in only about 0.5-0.8 additional photons 
being consumed in either simulation, in agreement with the lower 
overdensity of the LG region. Based on these results we therefore 
can conclude that in the photon-poor Model 2 the reionization of 
both the Local Group and of the nearby clusters is mostly a local 
process, with no significant external contributions from other struc- 
tures. 



4 SUMMARY AND CONCLUSIONS 

We have performed the first simulations of the reionization history 
of our local neighbourhood of the universe based on constrained 
N-body simulations of the formation of local structures - our Local 
Group of galaxies and the nearby galaxy clusters. The reionization 
history of the Local Group cannot yet be predicted uniquely, pri- 
marily due to the still poor observational constraints on the proper- 
ties of the reionization sources. Therefore, we studied two models 
constructed so as to roughly bracket the range of expected out- 
comes. While both models satisfy the available global observa- 
tional constraints, they differ significantly in their underlying as- 
sumptions. Our first model assumes relatively high ionizing photon 
production efficiencies ('photon-rich reionization'), while our sec- 
ond model studies the other extreme, where the sources produce 
barely enough photons to complete the reionization process in time 
Cphoton-poor reionization'). The second model therefore results 
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in more extended reionization history and delayed final overlap. 
The two models have the same observational constraints imposed 
on their initial matter distribution, which guarantees that the large- 
scale structures closely resemble our local neighbourhood, but they 
differ in their random component, yielding different constrained 
realizations. This allows us to check the robustness of our results 
with respect to the specific realization. Finally, the two models also 
have different underlying cosmological models (WMAP 3 -year vs. 
WMAP 5-year best fit). However, the effect of the background cos- 
mology is well understood and results in an overall shift of the 
reionization earlier or later, with no significant effects on our re- 
sults, for which only the relative timing of structure formation vs. 
reionization history is of importance. 

Our results show that the assumed efficiency of the ionizing 
sources has the most important influence on the nature of the reion- 
ization history of our Local Group of galaxies. Efficient photon pro- 
duction ensures that the nearby clusters emit more than sufficient 
number to ionize both themselves and their surroundings, includ- 
ing the Local Group. The fact that those galaxy clusters (Virgo and 
Fornax) coincide with high, rare peaks of the density field means 
that they form their progenitor halos earlier than the LG, which is 
in a more average region of the universe. As a result, the large-scale 
ionization fronts which propagated outward from the proto-clusters 
overrun the Local Group before it managed to form enough sources 
to ionize itself, resulting in its reionization being mostly externally- 
driven. 

Several points are worth noting here. Although generally the 
radiative transfer is a highly non-local phenomenon, which feature 
complicates its numerical treatment and the code parallelization, 
during most of the EoR the situation is somewhat more compli- 
cated. The neutral patches have enormous optical depth to soft ion- 
izing radiation (the only type of radiation we consider here). Even 
the already-ionized patches still have considerable continuum op- 
tical depth over cosmological (multiple Mpc) distances due to the 
small residual neutral fraction still remaining in such regiono This 
residual neutral fraction diminishes over time, but does so only 
gradually, as more and more sources appear and the mean flux 
thereby increases. As a consequence of all this, reionization starts 
out as a fairly local process where only the relatively nearby, di- 
rectly visible ionizing sources within the same ionized bubble con- 
tribute to the flux at a given point. This property allows us to focus 
our analysis on the important local sources and ignore the far-away 
ones for our current purposes (they are of course all included in the 
radiative transfer simulation). In our ^100 Mpc box there are mul- 
tiple proto-clusters which collapse nonlinearly by the present, but 
of those only Virgo and Fornax are sufficiently close to potentially 
contribute to the reionization of our Local Group. 

Furthermore, the ionization fronts propagate through under- 
dense regions (voids) much faster than through overdense ones (fil- 
aments, knots). Therefore, the relative positioning of the structures 
of interest and the density fluctuations in their immediate neigh- 
bourhood are important. Once the available observational con- 
straints are imposed in order to reproduce the local structures, we 
find that the Local Group is separated from Virgo and Fornax by 
voids in either realization (see Figure |2]). In contrast, the previous 
studies of this probleni which did not use c onstrained realizations 
JWeinmann et al .1120071 : lAlvarez et al.ll2009l) sampled a wide range 



^ E.g. At redshift z = 8, mean density and residual neutral fraction of 
10~^ the mean free path is ~ 3 Mpc. Density fluctuations further diminish 



this value. 



of environments and relative positions of nearby clusters. Such, 
purely statistical approach yields valuable insights on the range of 
reionization histories that could be expected for a certain type ob- 
ject (e.g. LG-like objects). However, by its nature such approach 
necessarily includes many objects which, although they share cer- 
tain basic features, locally do not reproduce the specific large-scale 
structures around us. Therefore, the constrained realizations are in- 
dispensible if we want to make realiable predictions for the effects 
of reionization on our neighbourhood. 

Why is the mode of reionization, external vs. internal, of our 
Local Group an important issue? This has a number of important 
implications for the formation of structures. Reionization dramati- 
cally rises the Jeans mass, thus impeding the formation and growth 
of small galaxies. In terms of this effect, the galactic haloes fall into 
three categories. The gas in the smallest halos (minihalos), whose 
virial temperatures are below the limit (^ 10^ K) for efficient ra- 
diative cooling through atomic line radiation. The ionization of the 
gas brings its temperature t o ^ 10^ K and it boils out, resulting i n 
their complete evaporation dShapiro et al.ll2004l : llliev et al.ll2005h . 
which leaves behind dark halos. In the other limit, the galaxies 
above certain mass (M > 10^° M©) have sufficiently deep gravi- 
tational potential wells to sucessfully withstand the effects of ion- 
izing radiation and are thus not significantly affected by the reion- 
ization process. The effects of radiative feedback on dwarf galax- 
ies of intermediate mass, roughly between lO^M© and lO^^M© 
is more complex and still very much a subject of active investiga- 
tion. The gas in such already-formed systems cannot be photoe- 
vaporated, as it can cool back down to ^ 10^ K very efficiently. 
However, photoionization heating rises the intergalactic gas tem- 
perature and pressure, which rises the Jeans mass and thereby sup- 
presses the future formation of very low-mass galaxies, as well as 
curtails the fresh gas infall onto such halos. Larger galaxies are less 
affected directly, but could do so indirectly, through their smaller 
progenitors, which could be expected e.g. to result in smoother gas 
sub- structure and modified stellar populations. Where the boundary 
between efficient and inefficient feedback from reionization lies is 
still unclear and very much subject of active research. Full investi- 
gation of the effects of reionization on galaxy formation and satel- 
lite galaxy populations goes well beyond the scope of the current 
work. However, our present results indicate that the photon produc- 
tion efficiencies of the first galaxies are the main factor determining 
the type of reionization history which our Local Group underwent. 
Therefore, this process should have left useful fossil records in the 
properties of our neighbourhood which will help us use local ob- 
servations to answer some of the key questions about the young 
universe. 
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